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Abstract 

■ 

We perform a joint jet tomographic analysis of the data on the nuclear modification 
factor Raa from PHENIX at RHIC and ALICE at LHC. The computations are performed 
accounting for radiative and collisional parton energy loss with running coupling constant. 
Our results show that the observed slow variation of Raa from RHIC to LHC indicates 
^ ' that the QCD coupling constant is suppressed in the quark-gluon plasma produced at 
LHC. 

1. The discovery of strong suppression of high-pT hadrons in AA-coUisions (usually called 
■ the jet quenching) is one of the main results from the RHIC program [1]. A similar effect 
^ ! has been observed in the ALICE experiment at LHC [2] for = 2.76 TeV. The most 
O ! natural physical reason for this effect is parton energy loss (radiative and collisional) in 
the hot quark-gluon plasma (QGP) produced in the initial stage of AA-coUisions. There 
has been much progress in the past 15 years in understanding the radiative energy loss 
due to induced gluon emission [3, 4, 5, 6, 7, 8]. For RHIC and LHC namely this mecha- 
nism dominates the energy loss [9]. Current studies in jet quenching are motivated by its 
importance for jet tomography of the dense QCD matter in AA-collisions. Unfortunately, 
^ I the available approaches to the induced gluon radiation are limited to one gluon emission. 

While for robust jet tomography, one should also take into account multiple gluon emis- 
sion. At present multiple gluon radiation is usually accounted for in the approximation of 
independent gluon emission [10]. However, this approximation has no serious theoretical 
justification in QCD, and can lead to considerable systematic errors. Also, additional 
uncertainties of the jet tomography come from the unsolved problem of treatment on an 
even footing the radiative and collisional (which is small but not negligible) contributions. 

Due to theoretical uncertainties presently the jet tomography can give only qualitative 
information about density of the QCD matter in AA-collisions at a given energy. How- 
ever, one may expect that the theoretical uncertainties should not be very important for 
variation of jet quenching with energy. For this reason the information extracted from jet 
tomography based on the data at very different energies should be more robust. From this 
point of view it is very interesting to perform a joint analysis of the nuclear modification 
factor Raa measured at RHIC and LHC. In the present work we perform such an analysis 
using the data on Raa for Au + Au collisions at y/s = 200 GeV from PHENIX [11] and 
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for Pb+Pb collisions at y/s — 2.76 TeV from ALICE [2]. A major purpose of this analysis 
is to decide whether the variation of Raa from RHIC to LHC indicates that the QCD 
coupling constant becomes smaller in the plasma produced at LHC. The suppression of 
the in-medium at LHC as compared to that at RHIC energies would be quite natural 
since the data on multiplicities [12, 13, 14] indicate that the initial entropy at ^/s — 2.76 
TeV is bigger by a factor of ~ 2.2 than that at ^/s — 200 GeV. This is translated into 
~ 30% growth of the initial temperature at y/s = 2.76 which should lead to a sizeable 
thermal suppression of a^. The fact that the PHENIX [11] and ALICE [2] data on Raa 
are very similar supports qualitatively suppression of at LHC. It would be interesting 
to see this from a quantitative analysis. 

The method of computation of the nuclear modification factor in the present work 
is similar to that used in [15]. We account for both the radiative and collisional energy 
losses. They are calculated with running frozen at small momenta. We treat the 
induced gluon radiation within the light-cone path integral (LCPI) approach [4]. The 
collisional energy loss is viewed as a perturbation effect. 

2. The nuclear modification factor Raa for a given impact parameter b can be written as 

^ dN{A + A^h + X)/dpTdy 

TAA{b)da{N + N + X)/dpTdy ' ^ ^ 

Here is the hadron transverse momentum, y is rapidity (we consider the central region 
y = 0), TAA{b) = J dpTA{p)TA{p ~ b), Ta is the nucleus profile function. The differential 
yield for high-p^- hadron production in AA-coUision can be written in the form 

dprdy J dprdy 

where dam.{N + N h + X)/dpTdy is the medium-modified cross section for the N + 
N ^ h + X process. Similarly to the ordinary pQCD formula, we write it as 

dam{N + N^h + X) _^ f^dz da{N + N^i + X) 

d^ -^Jo ^^'^A^^'^^ d^ ■ 

Here = pr/z is the parton transverse momentum, da{N + N ^ i + X) / dp\^dy is 
the hard cross section, D'^''^ is the medium-modified fragmentation function (FF) for 
transition of a parton i into the observed hadron h. For the parton virtuality scale Q 
we take the parton transverse momentum p^. We assume that hadronization occurs 
outside of the QGP. For jets with E 100 GeV the hadronization scale, jih, is relatively 
small. Indeed, one can easily show that the L-dependence of the parton virtuality reads 
Q'^{L) ~ max {Q / L , Qq) , where Qo ~ 1 ~ 2 CeV is some minimal nonperturbative scale. 
For RHIC and LHC, when tqgp ~ Ra {tqgp is the typical lifetime/size of the QGP 
phase, Ra is the nucleus radius), it gives fih ~ Qo (for E ^ 100 GeV). Then we can write 

/■I dz' 

D'^/iiz, Q) ^ / —Dn,Az/z\ Qo)D]}iiz', Qo, Q) , (4) 

J z Z 

where Dh/j{z,Qo) is the vacuum FF, and DJ'j-{z' ,Qo,Q) is the medium-modified FF for 
transition of the initial parton i with virtuality Q to a parton j with virtuality Qo- For 
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partons with E <^ 100 GeV the typical length scale dominating the energy loss in the 
DGLAP stage is relatively small ~ 0.3 — 1 fm [15]. This length is of the order of the 
formation time of the QGP tq ~ 0.5 fm. On the other hand, the induced radiation stage 
occurs at a larger length range / ~ tq -j- tqgp- For this reason to first approximation one 
may ignore the overlap of the DGLAP and induced radiation stages at all [15] . Then we 
can write 

/■I 

DT/i{z, Qo, Q)= ^DJ}f{z/z', E{)Df,f^^^{z\ Qo, Q) , (5) 

where Ei = Qz', DJjf is the induced radiation FF (it depends on the parton energy E, 
but not virtuahty), and D^^f-^^^ is the vacuum DGLAP FF. 

We have computed the DGLAP FFs with the help of the PYTHIA event generator [16] . 
Our method of calculation of the in-medium FF via the one gluon probability distribution 
is described in detail in [15], and need not to be repeated here. We just enumerate its 
basic aspects: 

• The multiple gluon emission is accounted for employing Landau's method as in [10]. 

• In calculating the q q FF the leakage of the probability to the unphysical region 
of /S.E > E is accounted for by renormalizing the FF. 

• The normalization of the FF for g ^ g transition, which does not conserve the 
number of gluons, is fixed from the momentum sum rule. 

• We also take into account the q ^ g FF which is usually ignored. Its normalization 
is fixed from the momentum conservation for q ^ q and q ^ g transitions. Thus 
for quarks our FFs satisfy the flavor and momentum conservation. 

We calculate the hard cross sections using the LO pQCD formula with the CTEQ6 
[17] parton distribution functions. To simulate the higher order i^-factor we take for the 
virtuahty scale in ag the value cQ with c = 0.265 as in the PYTHIA event generator 
[16]. We account for the nuclear modification of the parton densities (which leads to some 
small deviation of Raa from unity even without parton energy loss) with the help of the 
EKS98 correction [18]. For the vacuum FFs we use the KKP parametrization [19]. 

2. One gluon induced emission has been computed within the LCPl formalism [4]. The 
formulas convenient for numerical computations of the induced gluon spectrum are given 
in our paper [20], to which the interested reader is referred. We take rriq = 300 and 
nig — 400 Mev for the quark and gluon quasiparticle masses. These values were obtained 
in [21] from the analysis of the lattice data within the quasiparticle model for the relevant 
range of the plasma temperature T ~ (1 — 3)Tc. To fix the Debye mass we use the results 
of the lattice calculations for Nf = 2 [22] which give the ratio Hd/T slowly decreasing 
with T ifin/T ft; 3 at T ~ 1.5Te, hd/T ^ 2.4 at T ~ 4Tc). However, the results for Raa 
are not very sensitive to the Debye mass. 

We evaluate the gluon spectrum with running a^. The details of incorporating the run- 
ning coupling constant in the LCPl formalism are described in [9] . We use parametrization 
of as with freezing at some value at low momenta. For vacuum a reasonable choice 
is a^^ ~ 0.7. This value was previously obtained by fitting the low-x proton structure 
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function F2 within the dipolc BFKL equation [23]. A similar value of follows from 
the relation /q ^^^dQ ""^^ ^ 0.36 GeV obtained in [24] from the analysis of the heavy 
quark energy loss in vacuum. In vacuum the stopping of the growth of at low Q may 
be caused by the nonperturbative effects [24]. In the QGP thermal partons can give an 
additional suppression of ct^ at low momenta (Q ~ 2 — 3T). To study the role of the 
in-medium suppression of as we have performed the numerical computations for several 
smaller values of We are fully aware that this procedure, in which the in-medium 
suppression of the coupling is enforced on the average, so to speak, via modification of 
one parameter at any plasma temperature is very crude. We leave more accurate 
calculations with a temperature dependent parametrization of as for future work. 

The coUisional energy loss is small as compared to the radiative one but not negligible 
[9]. In the present work we treat the effect of the coUisional energy loss on the nuclear 
modification factor as a perturbation within the method suggested in [15]. It consists in 
renormalization of the initial temperature for the radiative contribution to the in-medium 
FFs using the following equation 

A£;,„d(ro) = i^Eradin) + AE,ol{To) , 

where AErad/coi is the radiative/coUisional energy loss, Tq is the real initial temperature 
of the QGP, and Tq is the renormalized temperature. As in [15] in calculating AErad/coi 
we take for the maximum energy loss half of the initial parton energy. The renormalized 
temperature is not very sensitive to the choice of the maximum energy loss in AErad/coi- 

We evaluate AE^oi within a modified Bjorken method [25] with accurate kinematics 
of the binary coUisions (the details can be found in [9] ) . We use the same infrared cutoffs 
and parametrization of the coupling constant for the radiative and coUisional energy loss, 
which is important for minimizing the theoretical uncertainties in the fraction of the 
coUisional contribution. 

We also included the effect of the energy gain due to gluon absorption from plasma 
by fast partons. It is done with the same prescription as for coUisional energy loss by 
renormalizing the plasma temperature. However, the effect of the energy gain on Raa is 
practically negligible [15], and can be safely neglected. 

4. We describe the QGP in the Bjorken model [26] which gives Tq to = T^r. We take 
To = 0.5 fm. Note that our numerical results show that Raa is rather insensitive to 
the precise vahie of tq. It is physically due to the fact that induced gluon emission, 
dominating parton energy loss, requires typically a finite formation time which exceeds 
considerably the value of Tq. For a given impact parameter the entropy density has been 
evaluated in the Glauber model using the Woods-Saxon nucleus density with parameters 
as in [2], and a'j^j^ = 42 mb for ^/s = 200 GeV and aj^j^ = 64 mb for ^ = 2.76 TeV. In 
calculating the entropy distribution in the impact parameter space we take it in the form 
dS/dy = C[adNpart/ dfj -|- (1 — a)dNcoii/ drj], where the Npart and Ncou terms correspond 
to the soft and hard mechanisms. We take a — 0.85 obtained in the hydrodynamical 
simulations of 74A-collisions in [27]. The normalization of the entropy has been fixed 
using the entropy/multiplicity ratio S/Nch = 7.67 from [28]. Note that for Au + Au 
collisions at ^/s = 200 GeV this prescription gives the total entropy smaller by ~ 20% 
than that used in [27]. To simplify numerical computations for each impact parameter b 
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we calculate the in-medium FFs for a uniform distribution of the initial temperature Tq 
which was obtained by averaging the realistic entropy distribution computed within the 
Glauber model. For the 0-5% centrality bin in Au + Au collisions at ^/s = 200 GeV we 
use dNch/dr] = 687 [12], and in Pb + Pb collisions at ^/s = 2.76 TeV dNch/dr] = 1601 [14]. 
These multiphcities give Tq a; 300 MeV for the central Au + Au collisions at ^/s = 200 
GeV and Tq ^ 400 MeV for Pb + Pb coUisions at ^/i = 2.76 TeV. 

The fast parton path length in the QGP, L, in the medium has been calculated ac- 
cording to the position of the hard reaction in the impact parameter plane. For L > tqqp 
we treat the medium as a mixture of plasma and hadron phases with relative fractions 
defined from the decrease of entropy density as s oc 1/r [26] . To take into account the fact 
that at times about 1 — 2 units of Ra the transverse expansion should lead to fast cooling 
of the hot QCD matter [26] wc also impose the condition L < L^ax- We performed the 
computations for L^nax — 8 and 10 fm. The difference between these two versions is small. 

5. In Fig. 1 the theoretical Raa for production in the 0-5% central Au + Au collisions 
at ^/s = 200 GeV computed with a^^ = 0.7, 0.6, and 0.5 for chemically equilibrium and 
purely gluonic plasmas is compared to the PHENIX data [11]. The results are presented 
for purely radiative mechanism and with inclusion of coUisional energy loss and radiative 
energy gain. As was said above, the effect of the radiative energy gain on Raa is practically 
negligible. The growth of Raa for gluons in Fig. 1 is due to the q ^ g transition which is 
usually neglected. However, it docs not affect strongly the total Raa since for ^/s = 200 
GeV the gluon contribution to the hard cross section is small at pt > 15 GeV. In Fig. 2 
we compare our results for = 0.7, 0.5, and 0.4 with the ALICE data [2] for charged 
hadrons in Pb + Pb coUisions at = 2.76 TeV. 

As can be seen from Figs. 1, 2, the collisional energy loss suppresses Raa only by 
about 15-25%. For the equilibrium plasma the data for ^/s = 200 GeV can be described 
with ^ 0.6 0.7. The data for ^ = 2.76 TeV agree better with af ^ 0.4 0.5. It 
provides evidence for the thermal suppression of ct^ at LHC due to higher temperature 
of the QGP. The in-medium suppression of as as the physical reason for a qualitative 
similarity of Raa at RHIC and LHC was also recently discussed in [29]. However, the 
effect has not been investigated quantitatively. 

Note that in our previous analysis [15] of the PHENIX data [11] the agreement with 
experiment was better for ~ 0.5. The difference with the present study is due to a 
somewhat larger total entropy of the QGP used in [15] . In [15] we used the entropy from 
[27]. As was said, the total entropy obtained using S/Nch ratio from [28] is smaller by 
~ 20% than that in [27]. However, this circumstance is not important from the point of 
view the conclusion on suppression of at LHC, since we use the same S/Nch ratio for 
RHIC and LHC. 

6. In summary, we have analyzed the data on the nuclear modification factor obtained 
at RHIC in the PHENIX experiment on Au + Au coUisions at = 200 GeV [11] and at 
LHC in the ALICE experiment on Pb + Pb collisions at ^ = 2.76 TeV [2]. The PHENIX 
data may be described in the scenario with chemically equilibrium plasma without (or 
small) in-medium suppression of ag- However, the data on Raa from ALICE, which are 
qualitatively similar to that from PHENIX, agree better with the results for ag smaller by 
~ 20 — 30% than that obtained from the PHENIX data. Thus, a relatively slow variation 
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of Raa from RHIC to LHC energies indicates that the QCD couphng constant becomes 
smaller in the hotter QGP at LHC. 
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Figure 1: The factor Raa for 7r° production in the 0-5% central Au + Au colhsions at 
^/s = 200 GeV for = 0.7, 0.6, and 0.5. The upper panels are for the chemically 
equilibrium plasma, and the lower ones for purely gluonic plasma. Black line: the total 
radiative part (quarks plus gluons); red line: the radiative quark part; green line: the 
radiative gluon part; blue line: the radiative (quarks and gluons) plus collisional, and 
plus energy gain due to gluon absorption. The theoretical curves obtained for Lmax = 8 
fm. The experimental points are the PHENIX data [11]. 
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Figure 2: The same as in Fig. 1 for the charged hadrons in Ph + Ph colhsions at y/s = 2.76 
TeV for = 0.7, 0.5 and 0.4. The experimental points are the ALICE data [2], as in [2] 
the boxes contain the systematic errors. 
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